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Abstract 

A road map to understand the relation between the onset of the superconducting state with the 
particular optimum heterogeneity in granular superconductors is to study a Random Tranverse 
Ising Model on complex networks with a scale-free degree distribution regularized by and expo- 
nential cutoff p{k) oc /c~'^ exp[— /c/^]. In this paper we characterize in detail the phase diagram of 
this model and its critical indices both on annealed and quenched networks. To uncover the phase 
diagram of the model we use the tools of heterogeneous mean-field calculations for the annealed 
networks and the most advanced techniques of quantum cavity methods for the quenched net- 
works. The phase diagram of the dynamical process depends on the temperature T, the coupling 
constant J and on the value of the branching ratio ^^^^^^^^^ where k is the degree of the nodes 
in the network. For fixed value of the coupling the critical temperature increases linearly with 
^^^^fc)^^^ which diverges with the increasing cutoff value ^ for value of the 7 exponent 7 < 3. This 
result suggests that the fractal disorder of the superconducting material can be responsible for an 
enhancement of the superconducting critical temperature. At low temperature and low couplings 
T and J ^ 1, instead, we observe a different behavior for annealed and quenched networks. In 
the annealed networks there is no phase transition at zero temperature while on quenched network 
we observe a Griffith phase dominated by extremely rare events and a phase transition at zero 
temperature. The Griffiths critical region, nevertheless, is decreasing in size with increasing value 
of the cutoff ^ of the degree distribution for values of the 7 exponents 7 < 3. 

PACS numbers: 64.60.aq, 64.70.Tg, 75.10Jm, 89.75.Hc 
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I. INTRODUCTION 



The interplay between disorder and superconductivity has attracted the interest of physics 
in the last decades. Disorder is expected to compete with superconductivity by enhancing 
the electrical resistivity of a system. In this situation by increasing the random disorder, the 
system undergoes a superconductor-insulator phase transition [ij. The theoretical explana- 
tion of the interplay between disorder and superconductivity is still a problem of intense 
debate for granular superconductors. 

Several authors have proposed that high-Tc superconductor are intrinsically inhomoge- 
neous {2 51 however these materials have a multiphase complexity that is difficult to tackle 
by analytical theoretical models in general. There is growing interest for a possible opti- 
mum inhomogeneity of superconducting cuprates that could enhance the superconducting 
critical temperature jo]. The control of defects and interstitials in heterostructures using 
new material science technologies can be used to design new granular superconductors with 
new functionalities js]. Recent experiments have provided multiscale imaging of the gran- 
ular structure of doped cuprate perovskites with a scale invariance that is reflected on the 
scale- free distribution of oxygen interstitials I?, 8|. Therefore it is now of high experimental 
interest the synthesis of novel granular superconductors made of superconducting networks 
with power law connectivity distribution. A road map to understand these new possible 
materials is to construct heterogeneous mean-field models and study the superconductor 
insulator phase transition in complex scale-free networks. Recently the superconductor- 
insulator phase transition has been characterized on Bethe lattices by solving the Random 
Transverse Ising Model on quenched Bethe lattices IJIllI with the use of advanced quan- 
tum cavity methods that allows to go behind the mean-field prediction. These methods 
have been recently developed to study exact phase diagram of the Random Transverse Ising 



model 
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13| . In this paper we use the quantum cavity mapping approximation of these 



methods, proposed in |9l-lll|. that gives a physical interpretation of the phase diagram of 
the Random Tranverse Ising Model by mapping the equation determining the phase tran- 
sition to a random polymer problem IJ]. In this paper we extend the results to networks 
with arbitrary degree distribution. In particular we focus of scale-free degree distributions 
mimicking the fractal disorder reported in cuprates. Already a mean-field calculation of the 
Random Transverse Ising Model on complex networks has shown that the superconducting 
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critical temperature is strongly enhanced on scale-free degree distribution with power-law 



exponent 7 < 3 



151 ] ■ This result might explain the experimental findings of Fratini et al. 



and might open a new road map to design new heterogeneous superconductors with en- 

t is in line 

with the study of other phase transitions _on scale-free network topologies 
was shown that the classical Ising model 
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spreading 



20|, the percolation transition 



17l |. In fact it 



2l| the epidemic 



22] change significantly when the second moment {k{k — 1)) of the degree distri- 



bution diverges. Recently it is becoming clear that these effects of the topology apply also 



to quantum critical phenomena [23 



d solution of both the 



24j . An open problem, 



and are found in the mean-fie 
Random Transverse Ising Model [l5j] and the Bose-Hubbard model 
that we approach in this paper is to what extent mean-field results are indicative of the 
behavior of the critical phenomena in quenched networks. In particular we focus here on the 
Random Transverse Ising model and we study the difference between the phase diagram on 
annealed and quenched scale free networks. While for annealed networks a simple mean-field 
results is guaranteed to be valid in quenched networks we have to approach the problem by 
the recently introduced quantum cavity method and by the mapping of this solution to the 
random polymer problem. We find that the phase diagram on annealed networks is signif- 
icantly different from the phase diagram of the quenched networks as long as the second 
moment of the degree distribution remains finite. In particular we found a phase transition 
at zero temperature not predicted by the mean-field approach and a replica symmetry broken 
phase a low temperatures. Instead, as the second moment of the degree distribution diverges 
{k{k — 1)) — 7- 00 the mean-field predictions approach the quenched solution. Therefore we 
identify the second moment of the degree distribution as a key quantity characterizing the 
disorder of the topology of the network. As the second moment of the degree distribution 
diverges, the superconducting critical temperature of the network diverges both on quenched 
and on annealed complex networks. 



II. RANDOM TRANVERSE ISING MODEL 



We consider a system of spin variables af, for i = 
given network with adjacency matrix a such that aij 



, . . . ,N, defined on the nodes of a 
1 if there is a link between node i 



n 



and node j otherwise aij = 0. The Random Traverse Ising Model is defined as in (SHuJ] and 
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le leading 
,0. The 



is given by the Hamiltonian 

ij i i 

This Hamiltonian is a simplification respect to the XY model Hamiltonian proposed by Ma 
and Lee {25! to describe the superconducting-insulator phase transition but to^ 
order the equation for the order parameter is the same as widely discussed in 
Hamiltonian describes the superconducting-insulator phase transition as a ferromagnetic 
spin 1/2 spin system in a tranverse field. We propose to use this Hamiltonian to describe 
in a granular superconductor in network with heterogeneity of the degree distribution. The 
spins (jj in Eq. ([T]) indicate occupied or unoccupied states by a Cooper pair or a localized 
pair; the parameter J indicates the couplings between neighboring spins, are quenched 
values of on-site energy and h is an external auxiliary field. To mimic the randomness of 
on-site energy we draw the variables ej from a p(e) = 1/2 distribution with a finite support 
e G (—1, 1). Finally in this model the superconducting phase corresponds to the existence 
of a spontaneous magnetization in the z direction. 

III. ANNEALED COMPLEX NETWORKS 

We consider networks of nodes i = 1,. . . ,N. We assign to each node i an hidden 
variable 6i from a p{6) distribution indicating the expect number of neighbors of a node. An 
annealed complex network is a network that is dynamically rewired and where the probability 
to have a link between node i and j is given by Pij 

p,, = P{a,, = 1) = ^ (2) 

In this ensemble the degree ki of a node i is a Poisson random variable with expected 
degree ki = 6i. Therefore we will have 



{9) = {k) 



{6') = {k{k-l)). (3) 

where (...) indicates the average over the nodes of the network and the overline in Eq. 
([3]) indicates the average over time- dependent degrees of the nodes. In order to mimic the 
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fractal background present in cuprates [8| we assume that the expected degree distribution 
of the network is given by 

p{e) = ue-^e-''^ (4) 

where is a normahzation constant and ^ can be modulated by an external parameter such 
as doping or strain in cuprates. 



IV. SOLUTION OF THE RTIM ON ANNEALED COMPLEX NETWORKS 
A. Critical temperature 

In oder to study the Random Trasverse Ising Model on annealed complex networks we 
consider the fully connected Hamiltonian given by 

ij i i 

where in order to account for the dynamics nature of the annealed graph we have substituted 
the adjacency matrix aij in H with the matrix pij given by Eq. (jj]). In order to evaluate the 
partition function we apply the Suzuki- Trotter decomposition in a number A^^ of Suzuki- 
Trotter slices. Therefore we have in the limit A^^ — oo 

Tre-^^ = Tr {e~^E/^^e-^^''^''^ ^^^^f' (6) 

where E is given by 

ij i a 

In order to perform this calculation we consider for each spin the sequence erf = 
{(jf'^, . . . crf'^°} where each spin af'" represents the spin i in the Suzuki- Trotter slice a. 
The partition function is then defined as 

N 

{2:Ji=i,...,iv *=i 
where we have indicated with 

a 
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In order to disentangle the quadratic terms we use A^^ Hubbard-Stratonovich transformations 

Ns/2 



\2tiNs ) J - 2N, 
exp[iV j d9p{9)^ j delnTr JJ 



where T>S_ = Yla=i (^S". The free energy / = hmTv-i^oo hm7Vs->oo j^^^^ can be evaluated 
at the stationary saddle point which is cyclically invariant. Therefore we get 

/ = inf^^^^-^ j d9p{9) c/eln(^2cosh(^/3v^(M~ZS^F^ 
where the value of 5* which minimizes the free energy is given by the saddle point equation 

« - E [ *7pi=^ ..MPVUseTWTI') (10) 

Finally the magnetizations along the axis x and z can be calculated by evaluating 

Tr a'e-^^ 



nrtn 



Z 

Tr afe-^" 



(11) 



Performing these calculations we get 

JS9 + h 



rrin 



rrin 



^{JS9 + hf + e2 



idJi\i{P^{JS9 + hf + e2 



^J{JS9 + h) 



2 _L £2 



Therefore the magnetizations m| ^ and m^' ^ depend on the value 9 of the expected degree of 
the node and on the on-site energy e. 

The order parameter for the superconducting-insulator phase transition is 5* given by Eq. 
ffTOj) . From the self-consistent equation determining the order parameter for the transition 
it is immediate to show that the superconducting-insulator phase transition occurs for /i = 
at 

(0) Jo e " MP) ^ ^ 

which implies that for {9"^) oo then /3 — j- for any fixed value of the coupling J > , and 

the critical temperature for the paramagnetic ferromagnetic phase transition Tc diverges. 
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This implies that on annealed scale free networks with p{6) oc O''^ and 7 < 3 the Random 
Tranverse Ising model is always in the superconducting phase. 

If we consider the special case of an expected degree distribution given by Eq. (jl]) and 
7 < 3 we have that diverges in the limit — ?■ 00. If we take ^ 1 which can be 

achieved by keeping J constant and we go in the limit — 00 with 7 < 3 we find that the 
critical temperature for the superconductor-insulator transition diverges as 

(9^) f Ine if7 = 3 
Tcoc J — — oc < 

{e) I ^3-7 if ^ < 3 

Therefore on complex topologies when the expected degree distribution is given by Eq. (j4]) 
and if 7 < 3 we observe an enhancement of the superconducting temperature with increasing 
values of the cutoff ^. On the contrary for Jj^r ^ 1 we have 



4e^ 

Tc = exp 

vr 



(0) 



(13) 



where C ~ 0.577 is the Euler number. Therefore in the annealed network model there is 
no phase transition at zero temperature. In fact from Eq. (fT3l) when (^^) is finite we have 
that Tc = only for J = 0. This phenomena strongly depends on the assumption that the 
network is annealed as it was found in M, uM where the critical behavior of the Random 



Tranverse Ising model was studied on a Bethe lattice. The critical indices of the phase 



transition will be given by the heterogeneous mean-field [19|] approach to phase transition 
as long as the cutoff = 00 and will depend non-trivially from the value of the power-law 
exponent of the expected degree distribution 7. The critical index of the susceptibility will 
be the only exception and remain mean-field as long as we don't include dependence of the 



network on the embedding space 20 1. 



V. QUENCHED NETWORK 

In this section we introduce quenched networks that have a structure that does not change 
in time. In the subsequent section we will characterize the Random Transverse Ising model 
on quenched networks and we compare the results with the heterogeneous mean-field results 
obtained in the previous sections. Therefore we consider a quenched network of nodes 
i = 1,2 . . . , N and degree distribution p{k). We indicate by N{i) the set of nodes neighbors 
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of node i. We assume that the network is random meaning that the probabihty to reach 
a node of degree k by following a random link is given by kp{k)/ (k). We assume that the 
distribution of quenched local energies ej is given by p(e) = 1/2 with e G (—1,1). We 
will consider in particular the case in which the degree distribution p{k) of the network is 
scale-free, with power-law exponent 7, and regularized by an exponential cutoff ^, i.e. we 
assume 

p{k) = Mk-'^e-^l^. (14) 
where A/" is a normalization constant. 



VI. SOLUTION ON A QUENCHED NETWORK 



In quenched networks critical phenomena might show a different phase diagram than 
on annealed networks. In particular this is true for the Random Transverse Ising model 
that on quenched Bethe lattices has a zero temperature phase transition which is absent 
on the corresponding annealed network as discussed in |9l-lll|. In order to study critical 
phenomena on quenched networks with a locally tree like structure (vanishing clustering 
coefficient) the theory of the cavity method has been developed and applied to a large variety 
of classical critical phenomena on networks. Recently this approach has been extended to 
quantum critical phenomena without disorder using the Suzuki- Trotter formalism finding 



exact results 
method 
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13| . In this paper we take a different approach by following a recent 



9l-lll| proposed to study the Random Transverse Ising model on Bethe lattices. 
This approach has the advantage that it can be applied to quantum critical phenomena 
with disorder and that it is analytically tractable. Moreover, as we will see, the solution of 
the problem can be mapped to a directed polymer problem shedding light on the nature of 
the different phases on the critical dynamics. 



A. Cavity mapping 

In the cavity model, the cavity graphs is considered. The cavity graph is the graph 
centred on a node i when one of its neighbor j G N{i) is removed. It is therefore assumed 
that all the other neighbors a G N{i) \j of node i are uncorrelated. On each of these nodes 
a it is assumed that an efficient cavity field 5^ j is acting such that the local Hamiltonian 
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FIG. 1: (Color online) Critical line Tc = Tc{J) for the Random Tranverse Ising Model on a complex 
network with power law exponent 7 = 2.5 and different values of the cutoff The dashed lines 
indicate the onset of the Griffith (replica-symmetry broken) phase. This picture show that the 

Grifhth phase occupies a reduced region of phase space if the cutoff ^ and the branching ratio 

(k(k-i) 
i,\ diverges. 




FIG. 2: (Color online) Critical value Jc{T = 0) for the insulator-superconductor phase transition 
and critical value J^^+^^^ for the replica-symmetric and replica symmetry broken phase. For 
7 = 2.5 the Grifhth phase occupies a reduced fraction of the phase space as the exponential cutoff 
^ diverges while for 7 = 3.5 the Griffith phase survives also in the ^ — > 00 limit. 
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FIG. 3: (Color online) The critical temperature of the superconductor insulator phase transition 
as a function of the cutoff ^ when the coupling constant J = 1 is kept constant. As 7 < 3 the 
critical temperature for superconductivity is strongly enhanced by the scale-free topology of the 
underlying network. 

acting on the cavity graph centred in node i is given by 

H^T=-'^^^- E M^ + ^-X + '/^.X]- (15) 

a€N{i)\j 

where we assume here and in the following that the external field h = 0. It was recently 
shown that the cavity-mean-field method is a very useful tool to study the Random Tranverse 
Ising model jsH]. In fact on the same time it provides an analytical approach to the 
solution of the Hamiltonian, and it provides a physical interpretation of the phenomena 
occurring at low temperatures. Finally this approximation is able to reproduce with good 
accuracy the phase diagram of the Random Transverse Ising Model. In the cavity mean-field 
approximation we replace the dynamical spin variables in H^'^" by their averages (a^) 
with values 

{<) = r^^^ tanh/3Ji?^ , + 4. (16) 

Therefore we assume that the cavity Hamiltonian acting on spin i in absence of spin j can 
be approximated by 

H-""^ = -e^at - Ja! (1^) 

a&N{i)\j 
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which imphes that Bij = J ^aeN{i)\j i'^a) therefore 



5. 



B,, = J I tanh/3y/i?g^, + 6^ (18) 



This recursion induces a self-consistent equation on the distribution Pk{B) of the cavity 
fields B for nodes of degree k 

HB-JY. - Jir^ tanh /3 Vi?,. + ] . (19) 
B. Mapping of the problem to a direct polymer 

In order to characterize the critical behavior of the Random Transverse Ising Model we 
can imagine to iterate the relation Eq. (fT8|) on the complex network for L ^ 1 times. For 
L finite and N ^ oo the corresponding graph is locally a tree with branching ratio 

The cavity field at the root of this tree is a function of the cavity fields on the boundary. 
If we assume that infinitesimal cavity fields 5 <^ 1 are acting on the boundary nodes the 
cavity field Bq at the root of the tree is given by 

^ = J2Il J^^^ = S (21) 

where V indicates the path on the tree from the root to the leaf of the tree, and the product 
riaGP is o'^^^ node along the path V. The critical point of the dynamics is given by 
the parameters for which In S = indicating when the cavity field propagates over large 
distances on the network. The expression Eq. ( 12T1) allows for a mapping of our problem to 
the directed polymer (DP) problem. In particular the function S is the partition function 
of the directed polymer where the energy on each edge is given by e-^" = and the 

temperature is set equal to one. This problem has been studied by Derrida and Spohn [l^ . 
In the directed polymer problem there are two regimes: 

• The replica symmetric (RS) regime in which the measure defined in Eq. (12T]) is more 
or less evenly distributed among the paths 
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• The rephca symmetry broken (RSB) regime in which the measure defined in Eq. ( 12T|) 
condenses on a small number of paths. 

In order to find the phase transition between the two phases we need to study the behavior 
of In S over a typical sample of the disorder. This is done by evaluating the average over 



the quenched variables of InH. This computation can be done by making use of the replica 
method 
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26| by computing 



InS 



lim ■ 

n-i>0 



n 



(22) 



where the replicated partition function is a sum over n paths 



E n(^ 



tanh /3ec 



(23) 



where the weight of each node a depends on the number of paths which go though this 
node. In the RS solution we assume that all the paths are non-overlapping and independent, 
therefore 



In J + In 



^ tanh/3e 
de 



(24) 



In the RSB solution we assume instead that the n paths consist of groups of n/x identical 
paths, and therefore = x. With this assumption we get 



—In S = In J H — In 

L X 



where 



/(^) 



In 



(fe(fe-l)) 
(k) 

{k{k-l)) 



de 
\ 



/ tanh /3e 



lnJ + /(x) 



X 



{k) 



^ , / tanh/3e 



(25) 



(26) 



The value of x is found by extremizing the function f{x) respect to x with x G [0, 1]. In 
particular by minimizing f{x) as a function of x one find that for /3 > /3_r5_b the solution is 
replica symmetry broken with x = m < 1 and for /3 > Prsb the function is minimized at 
the boundary m = 1 and the solution is replica symmetric. 



C. Phase diagram 

The critical behavior of the Random Tranverse Ising Model is found when -^In 2 = and 
when Bq/B the cavity field at the boundary propagates at long distances on the network 
under the iteraction Eq. (fTSjl . Therefore we have to distinguish between two regimes: 
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The rephca symmmetric (RS) regime for (3 < (3rsb- 



In this regime we have that the critical line is dictated by the relation jItlE 
ln(J) = 0, therefore the critical points are solution to the equation 

.(fc(fc-l)) 



J- 



{k) 



tanh/3e 

ae 

e 



In the limit J 



(fc(fc-i)) 
(fc) 



^ 1 we get 



(27) 



^ ^ (k) ^ 



In^ if 7 = 3 

^3-7 if 7 < 3 



where the last relation is valid if p{k) = Mk . 

The replica symmmetry broken (RSB) regime for /3 > Prsb- 



In this regime we have that the critical line is dictated by the relation -^In S = ln( J) + 
/(m) = and f'{m) = 0, therefore the critical points are solution to the system of 
equations 



r 



de 

V 



{k) 

f tanh /3e 



^ , / tanh (5e 

de 

V e 

m 



In 



J 



tanh fie 







(28) 



1. Phase diagram atT = 



It is interesting to study the critical point at T = 0. In particular we observe relevant 
deviations of the phase diagram on a quenched network respect to the behavior on a annealed 
network when the critical point is predicted to be at J = 0. The function /(m) defined in 
Eq. f l2B]) take the simple form 



f(m) = — In 
m 



{k{k-l)) 



(k) 1 — m 

At T = we can study the critical point of the transition by imposing f'{m) 



-j-lnS = ln(J) + f{m) = 0. By studying these equation we get the critical point 

{k{k-l)) ( {k{k-l)) 
Jc (k) V (k) 



(29) 
and 

(30) 



13 



Therefore on any given network with finite branching ratio there is a quantum phase transi- 
tion at T = 0. In the case in which p{k) = J\fk~"'e~''^^ and 7 < 3 we observe for large value 
of the exponential cutoff ^ that Jc goes to zeros as 



Jc oc 



if 7 < 3 
if 7 = 3 



In^lnQn^) 

Therefore in the ^ — 00 limit, as long as the power-law exponent 7 is 7 < 3 we recover the 
mean field result that there is no phase transition at zero temperature. 



2. Phase diagram for T > 0; RS-RSB Phase Transition 

It is interesting to find the critical point Tc = Trsb ^ 1 where we have the transition 
between the replica symmetric phase and the replica symmetry broken phase. This point is 
found by solving the equations Eqs. (128 p for m = 1 which we rewrite here for convenience. 



Jk{k-l)) , /tanh/3e 



/ tanh /3e 



In 



J 



e 

tanh f3e 



(31) 



We can find the critical point by expanding the equation for J-^^^j^y-^ ^ 1 and T <^ 1 
finding 



o (fc(fc-l)) i^ (^{k{k-l)) 



JRS-RSB \^ 



rpRS-RSB ^ ^ 



TT I {k) V (k) 

where C is the Euler constant C = 0.577... Therefore the critical temperature for the onset 
of the RSB phase for degree distributions given by Eq. (|T^ and 7 < 3 goes to zero as the 
exponential cutoff ^ diverges, and we have 

rpRSnSB ^ ) 52(3-7) in^(g) if 7 < 3 ^ 
[ [lnCln(lnO]^ if 7 = 3 

In Figure [1] we plot the critical temperature for the phase transition between the super- 
conducting and insulator phase for small value of the temperature T and coupling J. For 
temperature T < Trs-rsb the critical fine is dictated by the Replica-Symmetry broken 
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equations Eq. ( 13T]) . By plotting Tc = Tc{J) for networks with different value of the branch- 
ing ratio K = ^^"''^^^'^''^ we show that as K grows the replica symmetry broken phase shrinks. 
In order to show how severe it is this effect in Figure [2] we show the critical coupling constant 
jRSrSB ^YiQ critical coupling constant at T = Jc = Jc(T = 0) as a function of the 
cutoff and the power-law exponent 7. We show that as the power-law exponent 7 < 3 the 
branching ratio K diverges with diverging values of the cutoff ^ and in the limit ^ — )■ 00 we 
recover the mean- field results. Moreover as we show in Figure [3l as the coupling is fixed and 
J-^^^^y^ ^ critical temperature for the superconductor-insulator phase transition 

diverges with diverging values of the branching ration K = -^^^^^y^- Finally the ciritical 
indices of the phase transition will depend on the value of the exponential cutoff and the 
power-law exponent 7 of the degree distribution. Moreover, following similar arguments 



used in 



10| we cha show that the replica symmetric phase is a Griffith phase. 



VII. CONCLUSIONS 



In conclusion, we have studied the Random Transverse Ising model on networks with 
arbitrary degree distribution as a paradigm to study superconductor-insulator phase transi- 
tions. We have studied the model both on annealed and quenched networks using mean-field 
and quantum cavity method. In particular we have chosen a recently proposed approxima- 
tion of the quantum cavity method that allows for a full analytical treatment of the problem 
while characterizing the phase diagram going beyond the mean-field treatment. This method 
is based on a mapping between the cavity equation and the random polymer problem in 
quenched media. We have therefore characterized fully the differences between the phase 
diagram of the Random Transverse Ising model defined on annealead and quenched net- 
works. The model shows a significant dependence on the second moment of the degree 
distribution {k{k — 1)). In particular the superconducting critical temperature is enhanced 
in networks with greater second moment of the degree distribution. Moreover in the limit 
{k{k — 1)) — 00 the phase diagram of the quenched networks is well approximated by the 
phase diagram of the annealed networks. 
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